Quantitative evaluation of defect-models in superconducting phase qubits 
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We use high-precision spectroscopy and detailed theoretical modelling to determine the form of 
the coupling between a superconducting phase qubit and a two- level defect. Fitting the experimental 
data with our theoretical model allows us to determine all relevant system parameters. A strong 
qubit-defect coupling is observed, with a nearly vanishing longitudinal component. Using these 
estimates, we quantitatively compare several existing theoretical models for the microscopic origin 
of two- level defects. 



A key limiting factor of superconducting quantum co- 
herent devices is that they suffer from decoherence in- 
duced by their weak but non-negligible interaction with 
the environment! 1 -. The theoretical modelling of these 
interactions has greatly advanced our understanding of 
fundamental processes in the environment! 2 " and led to 
improved designs for increased coherence times, e.g., by 
engineering 'sweet-spots' or insensitivity to particular as- 
pects of the environment^. Despite these advances, not 
all effects of the environment are understood. One such 
enigma is the appearence of pronounced anticrossings in 
the spectra of superconducting phased and &\jo& qubits, 
which are indicative of a strong interaction with an ad- 
ditional quantum system. It has be shown that these are 
coherent 7 two-level, or at least strongly anharmonio^, de- 
fects, but their exact microscopic nature is still unclear. 

In several experiments^ — , it has been observed that, 
for strongly coupled defects, the coupling term is trans- 
verse (involving pure qubit-defect energy exchange) with 
minimal longitudinal (phase shift inducing) component. 
In this work, we perform a high precision comparison be- 
tween experimental data and a general theoretical model 
to shed light on the exact form of the coupling operator 
between qubit and two- level defect. We obtain quantita- 
tive estimates of the longitudinal and transverse compo- 
nents and then compare our results to existing theoretical 
models for intrinsic two- level systems. 

We theoretically describe the system of qubit and two- 
level system (TLS) by the Hamiltonian 

H = H q + H T ^ s + Hj (1) 

where H q describes the qubit, Htls the TLS and Hi 
the interaction between the two. Our qubit is a flux 
biased phase qubit£^, consisting of a superconducting 
ring interrupted by a Josephson junction and threaded 
by an external flux. The qubit Hamiltonian is given by 

H g = ^j-q 2 - Ej cos 4> + (^j ($- 0Ext) 2 , (2) 

where Ej = I c <$>o/2tt is the Josephson energy of the cir- 
cuit, C is the qubit's capacitance, L the inductance of 



the superconducting ring and $o is the superconducting 
flux quantum. Eq. ((2]) describes an anharmonic oscilla- 
tor with dynamical variables given by the phase differ- 
ence across the Josephson junction <p and its conjugate 
momentum q, corresponding to the number of cooper 
pairs tunneled across the junction, with [g, 0] = i. The 
external flux (j)Ext is generated via a flux coil on chip. 
We assume a linear flux-current relation of the form 
<\>Ext = c^bias + P, with the fabrication dependent pa- 
rameters a and j3. The TLS is described as a generic 
two-level system and we write its Hamiltonian in the 
eigenbasis i^TLS = ^tlsTs? with the level splitting 6tls 
and the Pauli- matrix r z . 

We consider three different coupling operators which 
may stem from fluctuations in the three terms of Eq. (|2]) , 
each of which corresponds to a different microscopic ori- 
gin. The state of the TLS may modulate the magnetic 
flux (j>Ext threading the superconducting loop 10 ! 11 or the 
critical current I c of the Josephson junction^ 2 -—, result- 
ing in coupling to (j) or cos 0, respectively. Alternatively, 
the TLS may couple to the electric field of the junction 
E oc g, which is consistent with the TLS being formed 
from a charge dipol o 15 i 16 . These three situations are de- 
scribed by the following coupling Hamiltonians Hi, 

H[^ = <f> (cos 0^ r x + sin #0 t z ) (3) 
Hj 2 ^ = v c cos 4> (cos 6 C r x + sin 6 C r z ) (4) 
H\ q) = v q q (cos q t x + sin q r z ) , (5) 

where v^, v c or v q parameterize the strength of the cou- 
pling. The angles 0</>,0 c ,0g G [0,7r] denote the relative 
orientation of the TLS eigenbasis, the physical meaning 
of which depends on the particular microscopic model. 

In order to compare the various coupling models, we 
define the transverse v±_ and longitudinal v\\ coupling in 
the qubit (|0), |1)) basis as 

2v ± = t; o cosfl o ((l|o|0) + (0|o|l)) (6) 
2v\\ = ^sin0 o ((l|o|l)-<O|o|O)), (7) 
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FIG. 1: (color online) Peak positions obtained for (a) qubit spectroscopy and (b) swap spectroscopy. For clarity, only 10% 
of the dataset is shown. The error bars give the 1-cr confidence interval for the fitted peak positions. The theoretical curves 
show the relevant transition frequencies for the coupled qubit-TLS system obtained via fitting the extracted peak positions (see 
text). Inserts show examples of the (normalized) escape probability as a function of excitation frequency and bias flux, from 
which the peak positions are extracted. 



where d = q, (f> or cos(/>, the qubit component of the 
coupling term given be Eqs. ([3j)- dSJ) . 

To shed light on the nature of the interaction between 
qubit and two- level defect, we need to determine the val- 
ues of v and 0. To this end, we have performed a series 
of spectroscopy experiments of a superconducting phase 
qubit strongly coupled to a TLS, at varying microwave 
power—, see Fig. [TJ Performing spectroscopy at both low- 
and high-power allows us to use a combination of single- 
and two-photon transitions to obtain spectral lines which 
are sensitive to the nature of the qubit-TLS coupling. 
We also performed 'swap-spectroscopy', where an addi- 
tional swap between qubit and TLS is performed before 
readout, effectively measuring the state of the TLS. We 
extract the frequencies of the various transitions in the 
coupled system by fitting each spectroscopic trace with 
Lorentzian functions. 

Our theoretical model, Eq. (pQ), can be described by a 
total of six independent parameters. Three parameters 
describe the qubit circuit and its tuning via the exter- 
nal flux: the critical current I c of the qubits Josephson 
junction and the parameters a and f3 describing the lo- 
cal generation of flux on chip and its coupling to the 
qubit loop. The TLS is described by its level splitting 
6tls and the interaction between qubit and TLS via v 
and 0. Figure [2] shows an illustration of the spectrum of 
the model and the influences of the different parameters. 
Since their effects on the spectrum, as indicated by ar- 



rows in Fig. [2j are all largely independent, this allows us 
to perform a fit to all six parameters simultaneously. For 
the circuit capacitance C and inductance L we take the 
design values of C = 850 £F and L = 720 pH. To account 
for fabrication variation, we repeated the fitting proce- 
dure with a ±5% tolerance in both L and C, resulting in 
no significant variation in the TLS parameter estimates 
(although 7 C , a and f3 vary accordingly). It is important 
to note that, since we are limited to spectroscopic data, 
our results are only sensitive to purely transversal oc g z t z 
and purely longitudinal oc cr x r x coupling terms. 

As an example, the estimated parameters for coupling 
to critical current according to Eq. (j4]) are: level split- 
ting e T LS = 7944.38 ± 0.08 MHz with coupling strengths 
v ± = 35.52±0.13 MHz and v\\ = 0.27±0.12 MHz (uncer- 
tainties correspond to 1-a confidence intervals through- 
out). We find the estimates obtained by fitting to each of 
the three coupling models are consistent with each other. 
Repeating the fitting for an additional defect in the same 
chip with different level splitting 6tls and coupling pa- 
rameters v, 6 produced qualitatively similar results, so 
we only consider one TLS in what follows. Full details 
can be found in the supplementary material^. 

We now discuss our results in light of several existing 
models describing the microscopic origin of such TLSs. 
Coupling to either magnetic flux or critical current gener- 
ates both a transverse and longitudinal component. Us- 
ing the ratio of these terms gives us an estimate for the 
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FIG. 2: Anatomy of a qubit-TLS anti-crossing in the high 
power regime. The overall slope of the spectral lines, their 
position and the spacing between one- and two-photon fea- 
tures allows us to calibrate the system, even for several in- 
dependent fitting parameters (see text). The separation and 
asymmetry of the lines, within the anti-crossing itself, allows 
us to estimate the transverse and longitudinal components of 
the coupling operator, respectively. 



orientation of tan# = 0.04 ± 0.02 for either coupling, 
placing strong constraints on critical current or magnetic 
flux coupling models. 

If the state of the TLS modulates the value of the 
magnetic flux threading the superconducting loop, the 
observed coupling would resuh— from a magnetic flux 
contribution of 5$ Ext = ^Ext^o/i^) ~ 250/i$ . As- 
suming the fluctuations result from a spin in the surface 
of the superconducting loop of wire-thickness ex l//ra, 
such a modulation of the magnetic flux corresponds to 
a magnetic moment approximately 10 5 times that of an 
electron spirJi 

In the model of Ku et al. 13 , a variation in critical cur- 
rent will couple to the qubit via the operator cos if). Using 
our estimates for v± and vy, we obtain^ a critical cur- 
rent variation of 81 c ~ 0.7 nA, where I c = 984 ± 2 nA. 
We can also use our estimates for tan# and 6tls to cal- 
culate the TLS Hamiltonian in its physical basis, where 
the two basis states correspond to different values of the 
qubits I c . We obtain Htls = 1/2 €qt z + 1/2 Aot x with 
e = 0.34 ± 0.16 GHz and A = 7.94 ± 0.01 GHz, giving 
two nearly degenerate states coupled by a large tunnel- 



ing element. A similar calculation also holds for magnetic 
impurities. 

Alternatively, the model of de Sousa et al^ assumes 
an impurity level in the junction which, via hybridisation 
with the Cooper-pairs in the superconductor, forms an 
Andreev bound state with energy inside the gap. Using 
this modeli 7 - results in an impurity energy of < 150 
MHz and a variation in critical current of 81 c < 1.5 nA. 
Such an impurity energy that is close to the Fermi edge is 
a consequence of the small longitudinal coupling, C w 0. 

For a purely transverse coupling to the electric field, 
following Ref. [H, we can estimate the (aligned) dipole 
size (as a fraction of junction thickness x) for our TLS as 
d/x = 0.08. Since the momentum operator q has no diag- 
onal component, this type of interaction would not lead 
to a longitudinal component ex a z r z in the coupling oper- 
ator. Spectroscopy therefore provides no direct measure 
of the orientation q of the charge-dipole. Determining 
additional components of this form requires experiments 
which probe the dynamical properties of the system^. 

Although the data is compatible with a small longi- 
tudinal coupling (fitting to flux or critical current cou- 
pling), the resulting coupling strength vu is comparable 
to the uncertainties and therefore we cannot rule out a 
pure charge-dipole. In such a case, a small longitudinal 
coupling component may also stem from a variation in 
the junction potential along the lines of Ref. 12, A linear 
combination of Eqs ((3j)- (J5J) is therefore also possible. 

Using general theoretical models and high resolution 
spectroscopy, we have estimated the various coupling pa- 
rameters between a superconducting phase qubit and a 
coherent two- level system within the qubit circuit. Com- 
paring with existing theoretical models, we obtained pa- 
rameter estimates for various suggested sources of such 
defects. In each case, the experimental data indicates 
a small or non-existent longitudinal coupling, relative to 
the transverse coupling term. These results allow us to 
place strong constraints on the parameters of the theo- 
retical models and test their validity. 
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